Conductivity in a disordered one-dimensional system of interacting fermions 
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Dynamical conductivity in a disordered one-dimensional model of interacting fermions is studied numerically 
at high temperatures and in the weak-interaction regime in order to find a signature of many-body localization 
and vanishing d.c. transport coefficients. On the contrary, we find in the regime of moderately strong local 
disorder that the d.c. conductivity ao scales linearly with the interaction strength while being exponentially 
dependent on the disorder. According to the behavior of the charge stiffness evaluated at the fixed number of 
particles, the absence of the many-body localization seems related to an increase of the effective localization 
length with the interaction. 
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The interplay of correlations and disorder in fermionic sys- 
tems is one of the challenging open questions in the solid 
state physics. The phenomenon of Anderson localization of 
single-electron eigenstates[ 1] is by now well understood in 
systems of noninteracting (NI) fermions. In particular, in one- 
dimensional (ID) systems all states become localized^ for 
arbitrary small disordery|] and hence there is no d.c. linear 
transport response at anytemperature T > 0. However, it had 
been long ago realized[4] that correlations among electrons 
as introduced via Coulomb electron-electron repulsion could 
qualitatively change transport properties of the system. 

So far, firm results and conclusions have been reached for 
the T — ground state of ID tight-binding fermionic system 
with a diagonal Anderson disorder. In particular, it has been 
shown by the density-matrix renormalization-group (DMRG) 
numerical studiesj^,)^] that in spite of correlations the many- 
body (MB) states remain localized, preventing the d.c. trans- 
port. The T > behavior appears to be much harder to 
deal with|7|, [8[] and, at present, the existence of MB local- 
ization beyond the ground state is controversial. ||9| The most 
interesting conjecture emerging from an involved analytical 
calculation! 10] predicts a finite-temperature phase transition 
between the MB insulator at T < T* and a conductor at 
T > T* . Since such a transition in fact implies a qual- 
itative change of character of MB states across the eigen- 
spectrum it as relevant and highly nontrivial to study sys- 
tems at high T — > oo . ill ill In this context, recent studies of 
energy-level statistics,! 1 1] the effective hopping in the config- 
uration space Il2ll and the decay of correlation functions lll3ll 
indicate a possible MB localization at very large disorder 
strength W.|12] The conclusions from the scaling analysis of 
the conductivity of such models appear similar. 1 14ll as well as 
the time evolution and the entanglement of wave-functions is 
concerned.! 15] 

On the other hand, recent direct numerical evaluation of the 
T > transport coefficients in disordered anisotropic XXZ 
model] 16] (model being equivalent in ID to a tight-binding 
fermionic system with nearest-neighbor interaction) does not 
show any indication of a crossover to a MB localization at 
low T or at larger W. This questions the conductor-insulator 



phase diagram and the relation to above mentioned studies. 

Our aim is to extend previous numerical studv lflrjll of trans- 
port properties of the ID disordered system, modeled by the 
t-V model of spinless fermions, in order to explore the phase 
diagram at high T with respect to the d.c. conductivity do- 
In contrast to most previous works, in which the interaction 
strength A = V/(2t) has been mainly kept fixed and pos- 
sible MB localization has been considered at large disorder 
values W, we start with a disordered system of NI electrons 
(A = 0), characterized by the vanishing d.c. transport at all T, 
i.e. (To = 0. By increasing gradually, at fixed W, the repulsive 
interaction A > we monitor a possible conductor-insulator 
transition in cr . Dealing with a finite-size system, instead 
of a singular behavior we expect that the insulator-transition 
transition should manifest itself as a crossover in oo vs. A. 
This crossover can be then used as a signature of a qualitative 
(gradual or abrupt) change of MB states with respect to the 
d.c. transport in the thermodynamic limit T — > oo. 

As the prototype model for the interplay of correlations and 
disorder we study the disordered ID t-V model. The Hamilto- 
nian represents a tight-binding band of spinless fermions on a 
chain, the repulsion occurs between nearest neighbors, while 
the disorder is in site energies, 

H = -i^(cJ +1 c J + h.c.) + V r ^n l+ in. i + ^e i n. i . (1) 

i i i 

By choosing site energies randomly in the interval — W < 
6i < W, we obtain in the NI limit V = the Anderson- 
localization model. In order to avoid the interaction-induced 
Mott-type insulator at A = V/(2t) > 1, we restrict our study 
to the regime A < 1 (note that for A = 1 the model (Q]l can be 
mapped on the isotropic Heisenberg model in a random field). 
We assume the chain with periodic boundary conditions and 
L sites. Furthermore, t = 1 is used as the unit of energy. 

To probe transport response we evaluate the dynamical con- 
ductivity cr(w), 
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with the current operator 

3 = itJ2( c i+i c i - h - c -) 



(3) 



We adopt the view that possible conductor-insulator transition 
needs to be a manifestation of the character of MB quantum 
state slfll , 12, f5] (hence not directly related to other thermo- 
dynamic quantities). Therefore one may as well restrict the 
study to the regime T — > oo, (3 — > 0, where all MB states con- 
tribute with equal weight to a(to). In this limit, the relevant 
(and nontrivial) quantity is a(uj) = Ta(uj). 

a(u>) is calculated by employing the microcanonical Lanc- 
zos method (MCLM) lll6ill7ll . particularly suited for dynami- 
cal quantities at elevated T. In the following we present re- 
sults for systems with L = 16 — 24 sites and for generic 
cases of half- and quarter-filling n = N e /L = 1/2,1/4, 
respectively. A sampling over N r ~ 100 random e% con- 
figurations is made to obtain the relevant average response. 
Finite size effects should not affect significantly our analy- 
sis in both, the energy and space domains. Concerning the 
space domain, we focus on disorder parameters W for which 
the Nl-electron localization length £o is much shorter than the 
size of the system, £o <C L. One may use the estimate (for 
V = 0) Co ~ 28.5/W 2 M which for W = 2 - 4 gives 
£ = 7 - 1.8 < L. Moreover, T = DMRG calculations in- 
dicate that the effective £ is further reduced by the interaction 
A > 0.H 

The energy resolution of our spectra is much smaller than 
the average level spacing associated with the largest system 
size L = 24 studied. Approximate eigenfunctions corre- 
sponding to T — > oo limit are converged in Mi ~ 2000 Lanc- 
zos steps, providing the energy resolution of SE ~ 0.004 (for 
L = 24). In the next step, AI2 ~ 4000 Lanczos iterations 
are used to evaluate a{bS), leading to an estimation of the fi- 
nal frequency resolution 5uj ~ 0.005. Since for the largest 
L = 24 the studied sector contains N st ~ 2 10 6 MB states, 
the average level spacing AE ~ 10~ 5 <C Scj, so the discrete- 
ness of the exact eigenspectrum due to finite system size plays 
practically no role in our results. 

In Fig. 1 we present typical high-T spectra for cr(cj), show- 
ing different W = 1.5, 2, 3, 4 for fixed A = 0.5. Since we 
deal with a substantial disorder, cr(w) are essentially differ- 
ent from the weak-scattering Drude-like form. All curves in 
Fig. 1 reveal maxima at u) m > moving up in frequency 
with increasing W, being a characteristic of a quasi-localized 
regime. At the same time, the optical sum rule is for T — > 00 
independent of W and A, [16] 



a(uj)duj = T^if) = nt 2 n(l - n), 



(4) 



whereas the d.c. value <7o = <t(0), being the central quantity 
studied further on, shows a pronounced variation with W. 

In Figs. 2a,b the emphasis is given to the low-w window, 
which is relevant for the extraction of the d.c. value do. 
W = 2 spectra for the n = 1/2 case are smoothed with 
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Figure 1. (Color online) Dynamical high-T conductivity a (to) for 
A = 0.5 and different disorders W — 1.5,2,3,4 evaluated for a 
half-filled system n — 1/2 and L = 24 sites. 
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Figure 2. (Color online) Low-oj part of a (to) for n = 1/2: a) A = 
0, 0.5, disorder W = 2 and different L = 16, 20, 24, b) W = 3, 
L = 24 and different A = - 0.5. 



an uj dependent damping 77 = rjo + {i]oc — ^o)tanh (lu/ojq), 
770 = 0.002 < Slu, r/oo = 0.02, io = 0.2. Such a damp- 
ing, used hereafter, preserves the sensitivity for lowest ui — >• 
frequencies, while smoothing configuration-dependent fluctu- 
ations (most pronounced at A — > 0) at higher ui. 

As argued above for £0 < L cases, Fig. 2a confirms the ab- 
sence of any evident L dependence (at least for A = 0.5). One 
can make an additional observation that, unlike for A = 0.5, 
the fluctuations of <r{iS) even at low frequencies remain sub- 
stantial for the NI case A = 0. One finds that these flue- 
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tuations diminish with the increase of the sampling N r over 
random disorder configurations, indicating that the repulsive 
interaction A > suppresses the sensitivity to the particular 
disorder configuration. In addition, A = 0.5 case in Fig. 2a 
reveals a remarkable linearity, ~ do + being ap- 
parently generic] 16] for all A > 0. 

By varying A = — 0.5, the role of the interaction on the 
MB localization is investigated in Fig. 2b. Results are given 
for fixed disorder W = 3, L — 24 and n — 1/2. Again, 
the remarkable linearity a(ui) at low-o; may be observed be- 
ing very reproducible for A > 0.1 in spite of very small d.c. 
values (To involved. Even more important, there is no signa- 
ture of a presumable qualitative change in a(uj) (at least for 
A > 0.1), which would point to the crossover for finite A 
from the MB localization, present at A = 0, to a conducting 
regime <t > 0. 

Based on the same numerical analysis of <t(w) used for 
Fig. 2a,b, in Fig. 3 we show the extracted d.c. values o-q vs. A. 
To suppress the effects of the configuration-fluctuating com- 
ponent of cr(w) as much as possible, ao is evaluated from a 
linear fit of <j(lo) in the frequency interval ui < 0.1. do ca l- 
culated in this way are given by symbols for W — 2,3,4, 
respectively. For A < 0.1, dashed lines are used in Fig. 3 to 
interpolate ctq to the theoretical value <5o = of the NI A — >• 
limit. 

Results in Fig. 3 are central to this work. In spite of small 
values of do, in particular for the W = 4 case, the extracted 
do show very consistent behavior. Namely, it is quite evident 
that in the interval W — 2 — 4 we do not find any signature 
of possible crossover in the behavior do vs. A, which could 
be interpreted as the onset of the MB localization for A < 
A C (W). In fact, the simplest dependence do oc A seems to 
represent well our results in the investigated regime 2 < W < 
4. 

Since the possible MB localization at A > should be 
more plausible at lower doping, where the condition £ < 1/n 
is stronger, we investigate the quarter-filling case n — 1/4 as 
well. Results for do vs. A are, however, even quantitatively 
similar to the n = 1/2 case, although, as expected, somewhat 
smaller values of do are obtained. As presented in Fig. 4, it 
is instructive to follow the dependence of the d.c. value do as 
function of the disorder W. We investigate the A = 0.5 case 
for two fillings, n = 1/2, 1/4, and the data for W = 1., 1.5 
are included. It is evident from Fig. 4 that the dependence 
is exponential, i.e. do ~ a exp(— bW), with b ks 1.7,2 for 
n = 1/2, 1/4, respectively. 

Fig. 4 gives also a clear limitation to our numerical ap- 
proach in the regime If > 1. Since the low-w slope a 
of d(ui) shows weaker dependence on W (a oc 1/W from 
Fig. 1), a reliable evaluation of do requires a very high res- 
olution (5u < 1. The latter is determined in our MCLM 
method by Mi , Mq, Lanczos steps, but in the final stage also 
by the MB level density and the average MB level spacing 
AE oc 1/N st oc exp(— (X). The reason for this is that macro- 
scopic results for the transport become plausibly relevant (for 
A > case) only if 5u> > AE. 
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Figure 3. (Color online) D.c. conductivity <jo vs. A for half-filling 
n = 1/2 and different W = 2,3, 4. Dashed lines for A < 0.1 are 
interpolations to the theoretically known value ctq (A = 0) = 0. 
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Figure 4. (Color online) ao vs. W for A = 0.5 and fillings n = 
1/2, 1/4. 



The above results show that a weak repulsive interaction in 
ID disordered tight-binding systems is capable of destroying 
the phenomena of MB localization. At the present stage of in- 
vestigations, we do not have a clear analytical or phenomeno- 
logical explanation for this property. However, we have so 
far discussed cases by keeping the fermion density n fixed. 
On the other hand, one can investigate the behavior of the MB 
states by fixing the number of fermions N e while changing the 
system size L. Several such studies have been reported. 11 811 
suggesting that the two-particle localization length £ in the 
presence of interaction becomes enhanced in comparison with 
the single particle case. Because the context investigated has 
been rather different, we study here this effect from the quan- 
tity directly relevant to the coherent charge transport, i.e. from 
the charge stiffness D. For T > 0, D is defined by[ 19] 



D = 



(5) 



As before, we are focused on the high-T regime /3 — > 0, when 
all the levels are probed, Z = N st . In the absence of disor- 
der and for finite A > 0, D — TD remains finite because 
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of the integrability of the model . lfl9ll With disorder switched 
on and N e fixed, one expects an exponential suppression of 
D with the increase of the system size L, D oc exp(— L/£). 
Furthermore, for the particular case of NI fermions £ should 
be independent of N e , i.e., £ = £q, with £o denoting the single 
fermion localization length. 

Since for the localized phase (e.g. for V > and T = 
0) the stiffness D is distributed for different realizations of 
the disorder according to the log-normal distribution, |5J we 
present in Fig. 5 the average D = exp((hxD)) vs. L, with 
N e = 2, A = 1 and various W = — 4. D is obtained by 
the full diagonalization (for N e = 2, N st oc I?). It is evident 
from the figure that, for larger disorders, W > 2, D decays 
exponentially with L, which is consistent with the MB local- 
ization in the n — > limit. On the other hand, it is remarkable 
that the interaction A = 1 increases D at large L. It seems 
from Fig. 5 that a crossover between the £ < £o an d £ > £o 
behaviors occurs for W ~ 2 and L* w 40. For strong disorder 
W = 3, 4, the crossover appears already at L* ~ 10. These 
results suggest that the increase of £ due to interaction may be 
an argument for the elimination of the MB localization at high 
T for finite fermion densities n = 1/4, 1/2 discussed in this 
work. 



not observe any clear sign of the latter, at least it is not evident 
enough. 

On the other hand, some recent numerical studies using 
different criteria, seem to point to the possible conductor- 
insulator transition and the MB localization at high T at much 
larger W, essentially within the same model with typically 
fixed interaction A ~ I. ll2l - ll4ll Translating the definitions 
of disorder W, their estimates for the onset of localization 
would be W > W* ^6 — 10, consistent with the observed 
qualitative change of the level statistics.] 11] It should be ob- 
served that such cases correspond to extreme disorder, which 
would require within our (or an analogous) approach the ob- 
servation (see Fig. 4) of ao < 10~ 4 . The corresponding res- 
olution Soj < 10~ 4 and large MB density-of-states leading to 
AE oc LW/N st < Sui may be in principle obtained e.g. by 
the full diagonalization for large enough L. However, the lat- 
ter is not reachable by up-to-date numerical methods. Hence, 
we cannot exclude such a scenario for the onset of the MB 
localization at W > W*, but on the other hand, such extreme 
disorder would also put limits to its theoretical as well as ex- 
perimental verification and relevance. 

We authors acknowledge helpful discussions with X. Zotos. 
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Figure 5. (Color online) Normalized charge stiffness D/N e vs. L 
for N e — 1,2 fermions with the interaction A = 1.0, for various 
W = - 4. 



In conclusion, our numerical results reveal a steady and uni- 
form increase of <r with the repulsive interaction A > 
at high T and at fixed disorder W, which does not sup- 
port a possible MB localization at A > in the considered 
regime. In order to put our results in broader context, let 
us comment the relation to other works. Authors advocat- 
ing a finite-T insulator-conductor transition [lOj] give an esti- 
mate for the transition temperature T* oc l/(7V£AlnA) where 
M is the (single-particle) density of states and A a character- 
istic matrix element for the electron-hole pair creation. For 
high T » 1 and weak interaction A < 1 one may translate 
the estimate \±M to a critical A* oc A* oc 1/(W£ ). Hence, 
for larger W > 2 we are in the regime where at least some 
crossover should be observed for A ~ A* > 0. Still, we do 
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